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■ The phenomenon of wave packet diffraction in space and time is described. It 
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consists in a diffraction pattern whose spatial location progresses with time. The 
pattern is produced by wave packet quantum scattering off an attractive or repulsive 
time independent potential. An analytical formula for the pattern at t ^ oo is derived 
both in one dimension and in three dimensions. The condition for the pattern to exist 



> 

. is developed. The phenomenon is shown numerically and analytically for the Dirac 

equation in one dimension also. An experiment for the verification of the phenomenon 
is described and simulated numerically. 
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1 Wave packet diffraction in space and time 

The diffraction process is a cornerstone of wave scattering. Quantum mechanical scatter- 
ing shows diffractive phenomena in space, hke nuclear diffractive scattering, and separately 
in time|[l|]. The combined effect of time dependent opening of slits for plane monochro- 
matic waves produces diffraction patterns in space and time.0 Both, diffraction in time, 
and diffraction in space and time are of the utmost importance in testing the validity of 
time dependent predictions of quantum mechanics. Recent measurements of atomic wave 
diffraction 0, have indeed demonstrated that the diffraction in time, process is supported 
by experiments. 

The present work will show that there exists a broad class of phenomena occurring in 
nonrelativistic and relativistic quantum wave packet potential scattering that behave as time 
dependent persistent diffraction patterns. The patterns are produced by a time independent 
potential. In the present case, the pattern will be shown not to decay exponentially, like for 
the diffraction in time phenomenon, and may provide a tool to deepen our understanding of 
time dependent quantum processes, pertaining to the description of atomic (or other) beams 
by means of wave packets. 

The peak structure exists for all packets, but, it survives only for packets that are initially 
narrower than a value related to the potential width. For wider packets, the peak structure 
merges into a single peak. 

The effect is named wave packet diffraction in space and time. It will be shown that the 
structure exists for repulsive and attractive potentials. It also appears when the relativisti- 
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cally invariant Dirac equation is used instead of the Schrodinger equation. 

The present work emerges from numerical investigations of one and two-dimensional wave 
packet scattering off an attractive potential that showed curious multiple peak structures 
resembling a diffraction pattern. Q It was seen there that, wave packets that are narrower 
than the well width initially, backscatter as a wave train that is coherent and multiple peaked, 
a diffraction pattern that travels backwards in one dimension and at large angles in two- 
dimensions. 

The effect is absent for packets whose initial width was much larger than the potential 
extension. For this case, a smooth wave hump proceeds both forwards and backwards. In 
section 2 we will treat the one dimensional case analytically. The Dirac scattering in one 
dimension will be dealt with in section 3. Section 4 generalizes the results of section 2 to 
three dimensions. In section 5 an experiment aimed at verifying the theoretical predictions 
will be suggested and simulated. Section 6 summarizes the paper. 

2 Wave packet diffraction in space and time in one dimension 

In order to develop an analytical formula for the difrraction pattern, we resort to the simplest 
possible case that can be dealt with almost completely analytically. This is the example of 
a wave packet scattering off a square well. 

Consider a gaussian wave packet impinging from the left, on a well located around the 
origin. 
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V{x) = -Vo e{w - \x\) (1) 

where, 2 wis the width, 9 is the Heaviside function, and Vq, the depth. The use of a square 
well facilitates the calculation, but, the results do not depend on the sharpness of the well as 
evidenced by the numerical results of 0]. (Theoretical arguments that support this statement 
may be found below after eq.(0).) Using the results of P, we can write immediately both 
the reflected and transmitted packets. From this point on we take h = c = 1. 

f°° _• _feL 

ipix,t) = / (f){k,x,t) a{k,q(i) e * 2 dk (2) 

J —OQ 

where x, t) is the stationary solution to the square well scattering problem for each k 
and a{k, qq) is the Fourier transform amplitude for the initial wave function with average 
momentum qg. 

Explicitly, in the backward direction that interests us here, for x < — w 



D{k, k') e' ^ + F{k, k') e-' ^ "'^ 



(^(/c, x, t) 
D{k,k') 
F{k,k') 
E{k,k') 

A{k,k') = {k + k'f e-'' ' ^' - {k - k'f e 



E{k, k') 
A{k,k') 

-2 i {k^ - k'^) sm(2 k'w) 

■J\2 -2 i y w 



l\2 „2 i k' w 



(3) 
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where 

a{k, go) = e-'^' (4) 

with k' = \Jk? + 2 m | Vq | , and a, the width parameter of the packet. 

The integral of the D term can be performed exphcitly. We will call this contribution to 

the wave ^pin■ For t —>■ oo we have 



TT 



a^ + it/{2 m) 



2 2 

777- , , o TTT- (J . I -9 /r-\ 

= ^77T(a;-a;o) -T— [X - xq- q^tlm) (5) 

For long times, the oscillations in the exponent of the term accompanying F in the wave 
function are extremely fast. The most important contributions come from phases that are 
an extremum with respect to k. We will call the contribution to the wave due to F, '?/'re/i- 
Using the saddle point method we find that, the extremal phase for '?/'r-e// demands the real 
part of the momenta to be positive for negative x, namely kextremai ~ 2 a2+f°7(2 m) 

We concentrate on the case k' >> k, the low energy regime, for which the polychotomous 
effect was most visible. [Q, ^ 

At /c ~ we have F ^ —1. To order j, for long times and distances x >> Xq » w 
using the properties of Gaussian integrals it is found that 



TT 

+ i t/{2 m) 



u = i {x + Xo)'^ - " {x + Xo + qo t/m)'^ (6) 



2t ' ' t 



The full wave in the backward direction is obtained from the sum of the interfering waves 
of eqs-d^J^). The amplitude of the wave for very long times becomes 



IV'I = 2 
= 2 
z = a' 

This expression represents a diffraction pattern that travels in time and persists. This 
result hinges upon the value of F as ^ 0. However, this value is essentially the reflection 
amplitude for zero momenta, that is equal to -1 regardless of the potential. Therefore, the 
expression for ip above does not depend on the choice of a square well. 

Figure 1 shows a comparison of the expression above (without any change in scale) to 
the numerical calculation of the reflected wave for t = 1.2 10^, a huge time compared to the 
transit time of the packet through the well. The initial momentum is go = 0.4, the well width 
w = 1, the initial position xq = —60, and the mass parameter m = 40. We chose a large 

mass in order to be on the nonrelativistic domain — = v = 0.01, for which the Schrodinger 

m 

equation is valid. The figure shows the theoretical expression and the numerical calculation 
agree, despite the approximate treatment of the F term. For the longest distances, the 
approximation is not valid. At these distances the approximation of A; ~ breaks down. 
Due to the same time dependence of both the sin and sinh pieces, it is quite clear that the 
pattern persists to infinite time. The larger xq, the cleaner the pattern, and more peaks 
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Figure 1: as a function distance x, numerical calculation (solid line) and theoretical 
expression of eq.(^), (dashed line) , at t=1.2 10''. The wave packet has average momentum 
q=0.4, mass m=40, width a = 0.5 and initial location xq=-60. The well parameters are: 
width w=l and depth Vo=l. 
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appear in the wave train. Clearly for xq — ^ oo, the pattern disappears. The smaller qo, the 
more visible -less background-, the diffraction pattern is. 

In order to understand the absence of diffractive structures for a broad packet we consider 
the imaginary term inside the sin function in eq.(|^. We have to estimate the extension of 
the diffraction pattern that changes with time. The momenta involved are kmax = — ^ "^"^^ 

t 

The constructive interference occurs whenever F of eq.(0) does not differ substantially from 
the value at k = 0. At higher momenta F is complex, the real part diminishes and the 
imaginary part starts increasing. This changes reduce the constructive part of the pattern, 
incoherent pieces arise that blurr the multiple peak behavior. The faster the drop of the real 
part of —F, the smaller the range of momenta for which the formula above holds. 

In figures 2 and 3 we have shown the real and imaginary parts of F for different well 
depths. Figure 2 corresponds to the case close to a transmission resonance at k = 0, 
Yj[ = 1.00658424209 and twice this value, with/^Q = ^J2 m |Vo|. Figure 3 corresponds 
to in between two transmission resonances at threshold ttSf— = 0.9998902939413 and one 
and a half times this value. Note however that the terminology of transmission resonance 
refers to a nonzero value of k. There is no transmission at zero momentum. In the actual 
packet scattering process, even if k'^ w is equal to the value of a transmission resonance or in 
between resonances, k' differs even by a small amount from the threshold value. Therefore 
our pictures display the expected behavior of F. Even if for a specific value of k the condition 
of a transmission resonance is met, this set is of measure zero as compared to the continuum 
of values of k' entering the calculation of the reflected wave. 
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Figure 2: Negative Real and imaginary parts of the coefEcient F of eq.^ as a function of k 
for the case ^ = 1.00658424209, w=l, and twice this value, w=2 
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Figure 3: Real and imaginary parts of the coefEcient -F of eq.ij^ as a function of k for the 
case = 0.9998902939413 w=l, and 1.5 times this value, w=1.5 
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We can see from the figures that in either case, close to a transmission resonance or in 
between resonances, the drop in —F scales as ^. The same is true for any other choice of 
the product /cq w. Therefore we can estimate kmax ~ (2 w)~^. With a proportionality factor 
that depends on the depth of the potential only. As it is our aim to understand the absence 
of a polychotomous wave train for a wide packet and not the exact cut-off, we proceed with 
the above estimate as a working hypothesis. For the imaginary part of the argument in 
the sin function in eq.(|^ to stay small and not blurr up completely the pattern we need 
2 kmax 0"^ Qo << 1, or, using our scaling argument 



For a fixed incoming momentum, the initial packet width has to be smaller than a certain 



interference between incoming and reflected waves D and F terms in eq.(^) is responsible 
for the diffractive pattern. The transmitted wave has no such two components. There is no 
possible interference and no diffraction peaks may be observed in the transmitted packet. 

All the above arguments pertain only to the specific example presently considered. How- 
ever, the example carries the essential ingredients for any type of packet. One only needs to 
replace the form factor in eq.(|^), by the corresponding expression for any other chosen 
packet. It does not depend on the type of potential. Other functional dependencies may 
give more complicated expressions than the simple sin function above, but essentially they 
will show the same structure. 




(8) 



proportion of the well width. This was indeed observed in numerical simulations [Q. The 



11 



The all important contribution to the diffraction in space and time comes from the 
interference between the incoming packet and the reflected wave originating from momenta 
k 0. The latter amounts to the excitation of a virtual state at almost zero energy quasi- 
bound state inside the well, as was indeed observed numerically in 0,0. The extreme 
importance of zero energy metastable states inside the well, is recognized in the literature 
for quite some time. In the limit of /c = 0, the virtual state becomes a half bound state. 
This state is termed usually a half-bound state due to a factor of one half in the phase shift 
at threshold in the derivation of Levinson's theorem in the nonrelativistic case[0, and the 
relativistic case of the Dirac equation.]^. Its existence affects the phase shift at threshold. 
The present work shows that the existence of metastable almost zero energy virtual states 
makes itself evident through wave packet diffraction in space and time. The effect should be 
independent of the sign of the potential. Transforming to a barrier is like transforming a slit 
to an obstacle (Babinet's principle). Both show diffraction in space. Therefore, we should 
see diffraction in space and time for barriers too, provided the width of the barrier is smaller 
than the initial width of the packet. This is indeed the case as shown in figures 4 and 5. 

A minimal uncertainty gaussian wave packet traveling from the left with an average speed 
V, initial location xq, mass m, wave number q = m v and width a, 

/ (x — Xq)^ \ 

i/j = C expii q {x - Xq) — j (9) 

is scattered from a Gaussian barrier 
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Figure 4: |^| as a function distance x for an initial wave packet of width a — 0.5 starting at 
Xo=-10 impinging upon a barrier of width parameter w=l and height Vq=0.2 at t=800, the 
initial average momentum of the packet is q—1. 
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Figure 5: |^| as a function distance x for an initial wave packet of width a — 2. starting at 
Xo=-10 impinging upon a barrier of width parameter w=l and height Vq=0.2 at t=1200, the 
initial average momentum of the packet is q—1. 
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2 

V{x) = Voexp(^-^^ (10) 

Let us now compare to the phenomenon of diffraction in tinie^ , and diffraction in space 

and time with plane waves0, The former arises when a shutter is suddenly opened at 

t=0. Behind the shutter there is a stream of monoenergetic particles. After the opening of 

the shutter, there arise oscillations in the current on the other side of it. These oscillations 

die out at large times. The latter is produced by the combined effect of shutters that open 

in time and interfere between them in space, or other variations[0]. Instead of having a 

vanishingly small time correlation length |1^, the process of wave packet diffraction in space 

and time, persists for infinite time. However, for the polychotomous pattern to survive we 

need an initial spread of the packet smaller than the extent of the potential. Translating 

to a framework of time scales, this would imply that the time scale of the opening of the 

shutter should be longer than the time scale generated from the width of the packet, namely 

Atshuttcr > — — , where, as above, a is the initial width of the packet. It appears interesting to 
Qo 

find out whether a finite opening time, may prevent the exponential decay of the diffraction 
in time process also. 

3 Diffraction in space and time for the Dirac equation 

The diffraction in time phenomenon |l[], does not exist in the relativistic case. This can be 
seen from eqs. (18) and (34) in ref.[|l|] for the ordinary wave equation and the Klein-Gordon 
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equation respectively. Only when the speed of light is taken to infinity, the nonrelativistic 
results are recovered, as can be seen from eq.(36) in the same work. This seems quite 
surprising because, the Klein-Gordon (KG) equation at low momenta compared to the mass, 
yields very similar phenomena as the corresponding Schrodinger (Sch) equation, as evidenced 
by a substitution ipxG = * i^sch- It is not possible to open a shutter in zero time, as 

it implies an infinite opening speed. A shutter that opens slow enough may provide a 
framework in which not only the pattern persists to long times, but allows its existence for 
the relativistic case in a causal manner too. 

In order to see that wave packet diffraction in space and time does exist in the relativistic 
case, we performed numerical calculations using the one- dimensional Dirac equation with a 
scalar potential Q S(x) as given in eq. ([IOD 



^ = 



For the initial packet we took a minimal uncertainty relativistically invariant wave packet 

^0 



u 

V 



dk e* ^(^~^'<')-2 (j2 {E Eo~k qo-m'^) 

i k 



dk 



k{x—X(,)—2 (E Eg—k qo—m^) 



E + m 



(12) 



^ The use of a potential in the non-stationary case for relativistic equations is questionable. However, it 
is still of relevance. 
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with E = y/k^~+ln?, Eq = ^ ql + m?, and we have taken 70 = , 71 = (Tx for the Dirac 
matrices, where a denotes the corresponding Pauh matrix. 



We extended the numerical method used for the nonrelativistic case[ll|, to the unitary 



evolution of the Dirac particle and found it extremely accurate. This method is a straight 
forward extension of the method of |Tl[|.^. The vector density p = J dx {\U\'^ + \V 



is 



conserved to an accuracy of more than 0.1% after millions of iterations. The equation is 
solved with high accuracy as determined not only by substitution, but, also by comparing 
to the free evolution of the wave packet above at large times. The price one has to pay for 
the Dirac case is a much smaller time increment because of the oscillations introduced by 
the mass parameter. Figures 6 and 7 show the results obtained. 

It is clear from the figures that the phenomenon exists in the Dirac case also. 

We now proceed to study the analytical behavior of the diffraction process for the Dirac 
equation. Consider the one-dimensional Dirac equation with a scalar potential only S(x) as 
given by eq.([TT|) For the initial packet we use a minimal uncertainty relativistically invariant 



wave packet of eq.(|T2|) The stationary wave function inside and outside the well for fixed 
momentum is given by (j){k, x, t) 



( 1 „-j E t \ 

1 e 



V 



/. j E t 

'2 e 



E + m 



numerical algorithm and routines are available from the author 
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Figure 6: Upper u—\U\ and ten times lower component v—\V\, of the Dirac wave function 
for an initial wave packet of width a = 0.5 starting at xq=-10, qQ=l on a well of width 
parameter w—1 and depth Vq—I at t—800 
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Figure 7: Upper u=|C/| and ten times lower component v=\V\, of the Dirac wave function for 
an initial wave packet of width a = 2 starting at xq=-10, qQ=l on a well of width parameter 
w^l and depth Vq^I at t^l200 
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bi{—w < X < w, k,t) 
b2{—w < X < w, k,t) 
4>i{w < X, k, t) 
4>2{w < X, k, t) 



i k' 



E + m* 

i k 
E + m 



{C e 



i k' X 



D e 



—i k' x\ 



F e' 



k X 



(13) 



with 



where ei = e' 



k w 



, 62 = e 

el (1 - gf - (1 + g)', k' = ^ E^ - (m*)2, m* = m - ^,.0 
The solution to the scattering of the packet becomes \l/(x, t) 
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—2 ei 64 1 


i + g) 


A 




2 ei 63 (1 


-9) 


A 




-Ag 




A 
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64 = 



k' w 



U{x,t) 
V{x,t) 



k' {E + m) 

k {E + m*) 

\ 
J 



(14) 



A 



U{x,t) 
V{x,t) 



(15) 



As in section 2 we consider the long time behavior of the wave. Here too, the reflected 
wave will receive contributions mainly from low values of k. Inspection of the reflection 



•^For a vector potential the * is in the energy instead of the mass 
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parameter B in eq.(|T4D, we find that for k ^ 0, B ^ — 1. The same liappened in 
tlie nonrelativistic treatment of section 2. For low initial momenta, go << rn,or v << 1, 
the exponential of eq.(|l|) becomes f» e" ™ 'h^^-''^ {qo-kf _ except for the 

factor *, that is independent of we are back at the nonrelativistic case for the upper 
component, and at the same time the lower component vanishes. Hence, the Dirac case in 
one dimension is exactly analogous to the nonrelativistic case for low average momenta of 
the initial wave packet. 

In the ultrarelativistic limit m ^ 0, the upper and lower components are identical up to 
a factor of i. The incoming and reflected waves do not spread. Hence, there is no interference 
between both. If we first find the wave in the backward direction as in section 2, and then 
take the limit of a small mass, we find that the pattern disappears, as the arguments of the 
sin function tend to zero. For t oo, the backwards scattered wave becomes a receding 
wave packet located around x = —Xq — t — w, without any diffraction. The incoming packet 
is no longer present in the region behind the well. 

In summary, the diffraction in space and time with wave packets exists in the relativis- 
tic regime, but gradually diminishes as the mass decreases, until is eventually washed out 
completely for the massless case, as was found for the diffraction in time process with the 
Klein-Gordon equation. ^ 
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4 Wave packet diffraction in three dimensions 

In section 2 it was shown analytically that there exist a broad class of phenomena occurring in 
nonrelativistic and relativistic quantum wave packet potential scattering that behave as time 
dependent persistent diffraction patterns. The patterns are produced by a time independent 
well or barrier. Differing from the diffraction in tinie^ and diffraction in space and time0 
phenomena, that apply to the time dependent opening of slits with plane monochromatic 
waves, the pattern does not decay exponentially with time. 

We now consider the generalization of the analytical formula of eq. (|^ to three dimensions. 
We take a gaussian wave packet in three dimensions, 



0(r, t) = 

yo = « qo ■ (r - ro) . (16) 

4 (T^ 

with qo, the initial average momentum of the packet, Fq, the initial location and a the 
packet width, impinging on a spherically symmetric square well 



V{x) = -Vo e{w - r) (17) 

where, w is the well width, and Vq, the depth, and G, the step function. The results 
are not dependent on the choice of wave packet however, the Gaussian profile facilitates the 
integrations. The total wave function with outgoing boundary condition at large distances, 
reads 
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^+ (F, t ) = ijj-^ (F, t) + ijjscatt (r , t) 

= J d^k /(k, qo) (^,„(k, r, t) + iPscattiK r, t)) (18) 



where 



k- r 



^l^scattik, r, t) = ^ ^ (2 L + 1) e^^^ szn{6L) P^r ■ k) 



i k r 



k r 



L 



/(k, qo) = e-e^-'^")' ^ * (19) 

The incoming wave may be found exphcitly to be 



= -4(zt/(fm)+a2)-'^°"' ^2°) 
Where D= |r — Fq — 2i cr^qol, erfc represents the complementary error function of 
complex argument and the absolute value in the definition of D pertains to the vectorial 
character, and not to the complex number. Error functions and similar Fresnel integrals are 
found in many diffraction formulae, such as those appearing in the works of Moshinsky. [jl[] 
Performing the angular integral over k, the scattered wave becomes 



i}scau{r,t) = J k^ dk $(r,/c,t) 
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$(f,fc,t) = 74^ ^ ^'^Hk-^'iD-i ^ (2 L + 1) e^^-.zn(5^)/i,M(A^,qo,?o,?)(21) 



Where 



Ll' ,m,m' 



J(2 / + 1) (2 /' + 1) 
= < / /' L|0 > <U' L\m m' M > (22) 



'(2L + 1) 

Here F;™ represents as usual the spherical harmonics, ji denotes the spherical Bessel 
function for complex argument, and < / /' L\... > represents the Clebsh-Gordan coefficient. 

The expression for is quite involved. However, for L = it simplifies considerably. 
Fortunately, this is the only contribution we need at long times, except for very specific 
cases for which the phase shift of higher partial waves dominates. Fot t oo the integral 
in eq.(PT|) oscillates wildly, except for values of k that give an extremal phase 0. It is easy 
to show that in this case, only extremely small values of k enter the integral. However, for 



very small values of k, reads 12 



'"^(^'^ = -(2 /-I)!! (2 / + !)!! ^^TTTT 

= ^ (23) 



where x = k' w = y/c^ + 2m|Vo| w, and ji, denotes the spherical Bessel function of 
order 1. Except for very specific cases for which there is a resonance at a higher partial wave, 
only the lowest partial wave matters. Higher partial waves are suppressed by factors of the 
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w ^ ' 



form — ^^^^^^^= . For very long times we can then approximate the scattered wave 
by taking only the L = term in eg . (pi]) . For this case II becomes 



^0 = 47r jo{k d) 

d = |ro + 2ior2qo| (24) 

where again the absolute value pertains to the vectorial character of the variable inside 
the vertical bars, and not to the complex number. Therefore, jo is really a complicated 
mixture of Bessel and Hankel functions. 

Inserting eqs.(E3|,p3|) in eg. (PT|) and recalling the rules for the integration of the error 



function for complex variables, we find the transparent result 



"o/ TT \l __^2^2fr + d .2 r-d .2 



^ d \it/{2m) +a^J \2r ^ v i; 2r ^ v 2; 



2^ {it /{2m) + a2) 



A. = . , ^'-'^ (25) 
2^ (zt/(2m) + (t2) 



With uq = w {1 — ) , K = ^2 m |Vo|, and d is defined in eg.( 

The incoming and scattered wave interfere in the total wave. Noting that the scattered 
wave depends on the absolute value of r, while the incoming wave depends on r itself, it is 
clear that there will be a different behavior at forward and backward angles. 

For long times and long distances we can safely approximate the complementary error 
functions with the value erfc~ 2. In order relate to the results of the one dimensional case, 
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we will bring here the expression for the full wave, in the case of a packet impinging from a 
distance tq at zero impact parameter, intial momentum go? on a well located at the origin. 
Using equations (p0| , |25|) we find 



) 
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m r 



(ro + 2 z qoa^) + ia) 



e" sin 



t 



e 



uq {2 d+ Uo) 
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,2 



(26) 



This is a diffraction pattern that travels in time as found for the one dimensional case. 
Again the pattern gets blurred and forms a single peak when the imaginary part of the 
argument of the sin function is large. As explained in section 2, this amounts to the condition 
o" >> There is here the additional blurring effect of the phase a. For narrow enough 

packets the pattern persists to infinity. 

The results are independent of the impact parameter provided the initial packet is located 
at a much longer distance along the direction of qo as compared to the transverse direction, as 
found numerically in ref.Q. Also the pattern does not depend strongly on the initial energy, 
as long as it stays small compared to the well depth or barrier height. The formula of eq.(^Bp 
may be used to find the angular distribution, the scattering and total cross sections, etc. 
Such expressions do not seem to exist in the literature for packets, even for the simplest 
square well case dealt with here. 

The present treatment took advantage of the properties of Gaussian packets and the 
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Figure 8: as a function of distance for various polar angles (f) and in the plane z—0 for 
three dimensional scattering. (See text) 
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simplicity of the square well. However, they are general. One has only to replace the 
Gaussian form factor by the appropriate one in case and the phase shift for the square well 
by the one corresponding to the specific choice of potential. The case of long range potentials 
deserves however some care. 

In order to visualize somehow the three dimensional results we present in figure 8 a 
hypothetical scattering event of a cold neutron packet of 1 Fermi width, with initial velocity 
of v=0.02 c impinging along the x-axis on a well located around the origin of width 10 Fermi 
and depth 40 MeV. The initial position of the neutrons is taken to be Xq = -20 Fermi, yo (the 
impact parameter) = 2 Fermi and the time is taken to be t = 5 10^^ in units of Fermi also. 
The numbers are taken just for the sake of exemplification and are unrealistic experimentally 
regarding the initial width of the packet. For a more viable experimental setup see section 5. 
We chose the potential range and width to give a large scattered wave contribution. Figure 
8 depicts the wave amphtude obtained by adding equations (|20|p5D as a function of distance 
in meters and we have multiplied the wave function by a factor of {it /{2m) + cr^)^, in order 
to have ordinates with values around |V^|=1. We have chosen to show graphs for the plane 
z=0. The pictures are similar for other planes. The darkened elements in the figure are 
due to fast oscillations of the wave. For this very advantadgeous choice of parameters the 
diffraction in space and time in three dimensions exists at all angles and not only in the 
backward direction. Perhaps this is expected for a very low speed packet that essentially 
spreads around a position close to the origin, while a secondary scattered wave interferes 
with it. Directionality is here weak, in contradistinction to the one dimensional case. The 
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Figure 9: as a function of distance for various polar angles 4> and in the plane z—0 for 
three dimensional scattering. Well depth increased by 5 % as compared to the previous 
figure 
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effect depends strongly on the input parametes. In figure 9 we modified the parameters of 
figure 8 only slightly by increasing the well depth by 5%. The effects are dramatic. Wild 
oscillations in small regions of space arise. This topic will be studied in a later work. 

5 Numerical simulation of a suggested experiment 

Present advances in bose traps permit the handling of atoms at very low temperature and 
velocities. At such energies the atoms move as if they were a wave packet provided the 
random agitation due to thermal effects is not as crucial. However, still the large amount of 
interactions between the particles and of them with the environment causes a decoherence 
of the wave in quite a short time. Instead of resorting to a cold bose gas, we will then focus 
on a drop of liquid Helium. Although the experimental details are beyond the expertise of 
the author, a hypothetical setup will be described in the hope that it is not too far away 
from reality. 

We will imagine a setup that may be dealt with approximately as a one-dimensional 
system, namely a drop of liquid Helium put on a surface a short distance apart from an 
impervious wall and far from any other boundary, inside a adiabatic container. Gravity is 
supposed to hold the drop in place vertically. The wall serves as a potential barrier and we 
are only interested in the horizontal spread of the drop due to the interaction with the wall. 

Consider a drop of liquid helium, of around 1 cm^ in volume lying near a corner of Cesium 
coated plates perpendicular to each other may be appropriate for this purpose. Cesium is 
needed in order to prevent wetting[|13|. The helium, plate and environment should be at a 
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temperature below the A point, although it would be interesting to see the effect of normal 
to superfluid component of the Helium on the wave packet behavior, as a tool to study 
decoherence. Put a large area detector far behind the drop and the plate. The number of 
particles measured as a function of time at a fixed position would then be determined by 
the absolute value of the wave function at that point. The particles are repelled by the wall 
and the packet spreads at the same time. 

Figure 10 depicts such a case. We solved the one dimensional Schrodinger equation for 
a packet as described in section 2 , with a barrier of very large strength. We assumed a 
total number of particles of = 5 * 10^^, obtained from the tabulated density of liquid 
hehum and a volume of 1 cm^. We used a plate (potential) of width w=l cm and strength 
Vq — AeV ^ 10^^, in units of ^ = 1. We took a detector size oi dx — 1 mm, and recorded the 
amplitude of the wave multiplied by dx in order to calculate the relative number of particles 
in the same distance and at varying time. The number of particles as a function of time 
-relative number multiplied by the total initial number in the drop for a properly normalized 
wave initially- is shown in the figure. 

It appears that, due to the large number of particles to be detected a simple weighing 
technique might be feasible, in this way we also avoid interaction of the packet with the 
detector. 

It is perhaps too optimistic to expect the proposed experiment will work so cleanly as in 
the simulation, due to all kinds of effects inside the drop that lead to decoherence, such as 
production of internal excitations like rotons, vortices, etc. However, some remainder of the 
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Figure 10: Number of Helium atoms counted at a distance of 5 cm behind a plate of thickness 
w=l cm as a function of time in seconds 
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effect might still show up in the counter. If it does, it will be a triumph for the quantum 
mechanical description of matter waves by means of wave packets of macroscopic size. 

6 Summary 

In summary, we have found in the present work that, the phenomenon of diffraction of wave 
packets in space and time, is determined by the spreading of the incoming waves and the 
subsequent interference between incoming and scattered waves as well as by the physical size 
of scatterers. It appears then, that the theoretical and experimental study of wave packet 
scattering as it evolves in time, gives us a new tool for the investigation of structure and 
dynamics of atoms and nuclei and a new possible testground for the predictions of quantum 
mechanics. 
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